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ABSTRACT 



For several hours on 2008 April 19 the T Tauri spectroscopic binary DQ Tau was observed to brighten, reaching a maximum detected 
flux of 468 mjy and likely making it (briefly) the brightest object at 3 mm in the Taurus star-forming region. We present the light curve 
of a rarely before observed millimeter flare originating in the region around a pre-main-sequence star, and the first from a classical 
T Tauri star. We discuss the properties and nature of the flaring behavior in the context of pulsed accretion flows (the current picture 
based largely on studies of this object's optically variable spectrum), as well as magnetospheric re-connection models (a separate 
theory that predicts millimeter flares for close binaries of high orbital eccentricity). We believe that the flare mechanism is linked to 
the binary orbit, and therefore periodic. DQ Tau makes a strong case for multi-wavelength follow-up studies, performed in parallel, 
of future flares to help determine whether magnetospheric and dynamical interactions in a proto-binary system are independent. 
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1. Introduction 

Observations of pre-main-sequence (PMS) stars at millimeter 
wavelengths are used to study thermal emission from gas and 
dust in their circumstellar disks where planets may be forming. 
At these wavelengths the disks are optically thin and the inte- 
grated flux is related directly to the total amount of cold, cir- 
cumstellar material. But on short timescales, we show how the 
dominating mechanism for millimeter emission can shift to high 
energy, short-period processes. We present 3 mm continuum ob- 
servations of the low-mass, classical T Tauri star (CTTS) DQ 
Tau, which was observed to brighten for many hours on 2008 
April 19, reaching a maximum detected flux of 468 mjy (com- 
pared to 17 mjy in the days both before and after the event). 

A stellar flare at these typically quiescent wavelengths, like 
a solar flare in the millimeter, is believed to be the superposi- 
tion of a gyro-synchrotron spectrum with that from synchrotron 
emission, the latter pea king in the sub-millimeter to far-infrared 
dKaufmann et alj[l986l) . The first light curve for a strong mil- 
limeter flare observed toward a y oung stellar object (YSO), 
GMR- A in Orion, was reported by iBower et alj d2003l) and at- 
tributed to magnetic activity. In protostellar environments, sce- 
narios where magnetic fields can interact, re-connect, and accel- 
erate electrons to high energies are thought to include: single-star 
stellar flares as we see in the Sun; interactions due to changes in 
the magnetic field lines between a star and its accretion disk; 
and, colliding magnetospheres and coronal flares between com- 
ponents of a multiple-star system. 

While longer-term radio variability is well documented 
in YSOs, and can also be the result of magnetic ac- 
tivity like starspots, few outbursts have been observed 
(Feig elson & Montmerld[l999l) . DQ Tau is only the fourth flare 



detected toward a PMS object, and the first from a CTTS. In 
addition to GMR-A, our observation follows the re-occuring 
flares from the we ak-line T Tauri sta r (WTTS) V 773 Tau A, 
studied in depth by iMassi et all d2002l 120061 120081) . and a sin- 
gle 1.3 cm flare from a second YSO in Orion named ORBS by 
iForbrich et all d2008l) . 

In Orion, follow-up study of GMR-A confirmed a photo- 
spheric origin of gyro-sync hrotron emission around a single 
magnetically-active WTTS dFuruva et al.ll2.003h , whereas obser- 
vations of the deeply embedde d ORBS could only limit the emit- 
ting region to a radius < 2 AU (Forbrich et al. 2008). In contrast, 
the periodic flaring from V773 Tau A is the first extra-solar ev- 
idence for interacting coronal flares, attributed in this case to 
collisions near periastron be tween extended h elmet streamers 
from the binary components (Massi et al. 2008); a phenomenon 
closely analogous to the situation for DQ Tau. With a higher or- 
bital eccentricity (e = 0.556, vs. 0.3 for V773 Tau A), the DQ 
Tau system geometry predicts an overlap of its two protostellar 
magnetospheres at each periastron pass age without the need for 
extended structure (Mathi etTet alJll997l) . 



2. The DQ Tau System 

DQ Tau (a = 04:46 :53.06, S = +17:00:00.1) was found by 
Mathieu et a D (11997b to be a non-eclipsing, double-lined spec- 
troscopic binary, comprised of two relatively equal-mass stars 
(~0.65M Q ) with spectral types in the range of K7 to Ml and 
a robust orbital period of 15.804 days. At closest approach, the 
binary separation is 8R* (= 13 R ), roughly four times smaller 
than the 30 R* (=6OR ) separation for V773 Tau A, and less 
than tw ice the expected T Tauri stellar mag netosphere size of 
~5R, JShuetalJll994t iHartmann et aT]|1994l) . As a result of 
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the binary geometry and magnetospheric m odels, re-connectio n 
events have been predicted for this source (Math ieu etalJI 1997b . 
but never tested observationally. We provide evidence that this 
effect is occurring and that it may be repeatable and relevant to 
many similar T Tauri systems. 

The spectral energy distribution (SEP) of DQ Tau is rather 
typical for a CTTS. It is best fit bv lMafhieu et all d 1997b using 
a large circumbinary disk of about 0.002-0.02 M . They find an 
additional 5 x 10 -10 M of warm (~1000K) circumstellar mate- 
rial close to the stars is necessary to explain the near-infrared 
excess, but not enough to maintain the obse rved accretion rate 
of 5 x 10 -8 M yr -1 (see lHartigan et al.ll 1995b . DQ Tau must be 
transporting material into its inner regions. Both a broad bright 
10/im silicate feature (see lFurlan et al.l2 006). and the absence of 
strongly blocked forbidden line transitions in the receding stel- 
lar flow, support op tically thin material in a region up to 0.4 AU 
around the binary dBasri et al.lll997t iNaiita et alj|2003b . 

At optical wavelengths the system is known to brighten by 
0.5 magnitudes coinciding with periastron, either shortly before 
or at closest approach. The mechanism regulating the optical ac- 
tivity is variable, as brightenings have been observed for only 
65% of periastron encounters, and in iso lated events, multipl e 
brightenings were seen in a single orbit (Mathieu et "ail[l997b . 
All brightenings were recorded within an orbital phase window 
of 0.30 (or 5 days) centered around a phase (O) o f 0.9, where pe- 
riastron defines cD = 0.0, 1.0 dMathieu et alj 1997b . This evidence 
for accretion presents an additional puzzle, indeed for most bina- 
ries, since a large inner gap is expected to be cleared by the res- 
onant and tidal forces generated by stellar movements, thereby 
quickly isolating a circumbinary disk from one or two circum- 
stellar disks. The small binary separation of DQ Tau, however, 
reduces the likelihood of stable circumstellar disks. These effects 
combined should speed up disk dissipation and halt accretion 
processes. Instead, it is unclear how easily material can cross the 
predicted gap and why some binar ies show inner cle ared regions 
and others, like DQ Tau, do not (iJensen et al.ll2Q07b . 

In this context, DQ Tau was t he first system to b e studied in 
terms of pulsed accretion flows (Mathi eu et al.lll99~7h . Extensive 
observations of the system's optical variability and photospheric 
line properties show increased veiling effects, spectral bluen- 
ings, emergence of outflow signatures in spectral lines, and the 
double-line structure in Ca II suggesting that bright blue-shifted 
materi al is being emitte d near the surface of the star during an 
event (Ba sri et al]|l997b . All these effects at periastron are con- 
sistent with pulsed accretion fl ows based on the theoretical mod- 
els of Artvm owicz & Lu bow (1996) that describe how the ec- 
centric binary orbit periodically perturbs the outer disk, causing 
material to collect at stable points near the inner rim, and then to 
stream across the gap. However, the theory does not require the 
peak accretion rate to occur near periastron. Instead it has been 
postulated that the accretion mec hanism is further reg ulated by 
strongly variable magnetic fields (Mathi eu et al.lll99~7l) . 

3. CARMA Observations 

We observed DQ Tau at 1 15 GHz (« 2.7 mm) with the Combined 
Array for Research in Millimeter-wave Astronomy (CARMA) 
during 5 tracks in April 2008, totaling 4 hours on-source. The 
millimeter array, located in eastern California (USA), provides 
150 baselines of 30-350 m in its C configuration and records 
one linear polarization. Our program cycled through 5 Taurus 
sources during each track, scanning each object for 12 min- 
utes and obtaining 1-6 scans per source, depending on the track 
length. A 5-minute scan of the gain calibrator 3C111 was in- 



Table 1. DQ Tau Flux and Intensity Data at 2.7 mm (1 15 GHz) 



JD (2454500+^ 


Intensity (mjy beam ') 


Flux (mJy^_ 


76.3330 


37.4 + 


15.f£ 




76.3923 


345.0: 


±19.8 


419 


76.4534 


355.1: 


±11.4 


468 


76.5132 


216.8 


±9.2 


284 


76.5789 


150.3 


±6.1 


194 


Quiescent 


13.3: 


±1.5 


17 



" Indicates the start of the 12-minute scan. 

b Determined by a gaussian fit to the continuum image. 

' Data point suffers from severe de-correlation at longer baselines. 



eluded every 12-24 minutes. Multiple observing scripts allowed 
rotation of the source order in each track so that we could achieve 
more efficient (u, v)-coverage per source. In addition to the phase 
calibration using 3C1 1 1, we observed the planet Uranus for flux 
calibration at the start of each track, as well as the radio source 
3C84 for passband calibration. 

All data were processed using the MIRIAD data reduc- 
tion software program, optimized for CARMA. We calibrated 
each track separately before combining the data into a final 
(w, v)-averaged continuum image. In the inversion step a natural 
weighting for the (u, v)-coverage was applied, and cleaning was 
performed down to a cutoff of 2cr. The continuum maps are a 
composite of two 500 MHz bands covering the frequency ranges 
110.94-1 11.41 GHz and 1 14.80-1 15.27 GHz, corresponding to 
the lower and upper sidebands, respectively. 



4. Results and Analysis 

In the combined and calibrated DQ Tau image (~3 hours on- 
source) we achieved a sky rms noise level of 1.45 mjy beam -1 
and a beam size of 1.92" x 0.99" toward the Taurus star-forming 
region, corresponding to a spatial resolution of -140-270 AU at 
a distance of 140 pc. Since we found that the source was con- 
siderably brighter on 19 April than during any other observation 
date (Fig. [TJ, data for that track were not included in this final 
continuum image. And thus the quiescent values for the other 
nights combined then give a peak intensity of 13.3 mjy beam -1 
and a flux, determined by a gaussian fit, of 17 mjy. 

Review of the unexpected result from 19 April revealed 
that the brightening was a real effect. Many instrumental and 
data processing errors could be ruled out because none of the 
4 other sources exhibited a similar effect, the brightening per- 
sisted throughout 5 source cycles in an 8-hour period, and, the 
same observing script was used for 3 of the 5 tracks. The loca- 
tion of the peak emission at the time of the flare is consistent 
with the coordinates of the binary on all nights and therefore no 
background source is suspected. Finally, while the linear feeds 
do rotate while tracking the source, we expect a flare mimicked 
by the presence of strong polarization would not only be present 
in the nearly identical tracks of the adjacent days, but that the 
emission profile would be more symmetric about its peak. 

To understand the nature of the DQ Tau brightening, we re- 
duced each of the track's five 12-minute source scans separately 
to produce individual continuum maps. We then plotted the max- 
imum source intensity (Jy beam -1 of the peak pixel) as a function 
of time (Fig.[TJ. DQ Tau clearly exhibited a flaring behavior, ris- 
ing to a maximum brightness in ~ 1-2 hours, coinciding with the 
start of our track, and shortly thereafter, decaying (roughly expo- 
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Fig. 1. A plot of the maximum intensity of each scan as a function of Julian Day. Filled squares represent DQ Tau, circles the 2 
brightest Taurus sources in our tracks, and crosses are the calibrator (divided by a factor 50). Error bars for DQ Tau are the sky rms 
for each scan. Upper Left Inset: Continuum images of DQ Tau at 2.7 mm (115 GHz). White crosses indicate the 2MASS infrared 
position. Pixels are 0.2" x 0.2" and the same grayscale is used for both images (ranging -3.14 x 10~ 2 to 3.55 x 10~' JybearrT 1 ). On 
the left, the combined image from the 4 quiescent tracks. Contours are drawn at 3<x levels where cr- 1.45 x 10~ 3 . At right, DQ Tau 
brightens by a factor of 27 on 2008 April 19 (JD = 2454576.4534), peaking at 355 mJy beam -1 . We determine a flux of 468 mJy 
for this scan. Upper Right Inset: The intensity (Jy beam 1 ) of the peak pixel at the location of DQ Tau throughout the track on 19 
April, now sub-divided into 2-minute intervals. The horizontal dashed line indicates the maximum intensity found for the quiescent 
image. The second and third scans show a sharply increasing and decreasing intensity, respectively. 



nentially) throughout the remainder of the observation, returning 
to half the observed peak intensity within a 2-hour period. 

In Fig. Q] we show the results of performing the same in- 
tensity analysis on the individual scans for three additional 
sources (our two brightest Taurus sources plus the gain cali- 
brator). These sources appear constant in flux for all observa- 
tions (distributed over 13 days). The same is not true for DQ 
Tau, which remains unresolved throughout its brightening event. 
At peak brightness its flux during our third scan on 19 April 
(JD = 2454576.4534) was 468 mJy, roughly 27 times the quies- 
cent value (Table Q]), likely making DQ Tau the brightest ob- 
ject at 3 mm in Taurus at the time of our observation if we con- 
sider t hat the class-defin ing source T Tau is just 56 + lOmJy at 
3 mm dOhashi et al.lll996|) and the brig ht, extended source HL 
Tau is 100 mJy ( Beckwith et al. 1986). Finally, the flare is in 
stark contrast to (single-dish) measurement s of DQ Tau that find 
125 mJy at 1 mm and 91 m Jy at 1.3 mm dBeckwith & Sa rgent 
Il99lt iBeckwith et alJI 1990h . Our quiescent value of 17mJy, on 
the other hand, is a better fit to the expected millimeter spectral 
slope. 

The true shape of the flare onset and its exact peak are less 
precise. Larger error bars for this date are a result of significant 
flagging due to strong decorrelation on longer baselines, particu- 
larly during the first half of the track. The gain calibrator, which 
has an average flux of 3.7 Jy for each track, has a larger stan- 
dard deviation (cr) for the maximum intensity of each 5-minute 
calibrator scan on the flare date (0.257 Jy birr 1 ) than, for exam- 
ple, the longest tracks on 13 April (0.095 JybrrT 1 ) and 21 April 
(0.103 JybtrT 1 ). These effects do not compromise the validity 
of the flare detection. For comparison, the two brightest Taurus 
sources observed have a cr of O.OOJybrrT 1 and 0.006 JybrrT 1 
for their five 1 2-minute scans during the flare period. 



Our flare occurs within the orbital phase window of the 
known optical brightenings, although it is unknown whether DQ 
Tau flared simultaneously in the optical. Using a JDo value of 
2449582.54 ± 0.05 and a period of 15 . 8043 + 0.0024 days, as 
published in Table 2 of Mat hieu et al.l dl997l) . we calculate the 
phase of the flare peak to be 0.98, or ~7.6 hours before closest 
approach. If we consider the error in the orbital period, propa- 
gated over 1 1 years, the flare in fact could have peaked anywhere 
between 25.4 hours before periastron (O = 0.93) to 1 1 .0 hours af- 
ter (<I> = 0.03), still well within the window of optical variability. 
Our other observations began at orbital phases of -0.40, -0.09, 
0.11, and 0.42 with respect to the same periastron encounter 
(corresponding to -6.32, -1.42, 1.74, and 6.64 days to closest 
approach), but only the track nearest periastron (O = -0.02 on 
19 April) exhibits a brightening event. Therefore, we strongly 
suspect that the millimeter flare is connected to the binary mo- 
tions, as well as to the optical variability, suggesting that this 
event is likely to be repeatable, and even periodic. 



5. Discussion 

Of the three previous examples of flaring activity, DQ Tau draws 
the strongest parallel to V773 Tau A both in system geome- 
try and flare properties. The similarities include: binarity with 
a fairly equal mass distribution, high orbital eccentricity, and 
now an observed 3 mm flare that is nearly identical in rise time 
(~2 hours), decay time (~5 hours, if we extrapolate in the case 
of DQ Tau), peak intensity (~400mJy), incidence at periastron, 
and ratio (~25) of maximum intensity reached with respect to 
quiescent levels (see lMassi et alj|2006l Fig. 2). The most signifi- 
cant difference is the lack of ongoing accretion around V773 Tau 
A as determined by its Ha equivalent width. 
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In the other examples, the coronal flare observed toward 
GMR-A produced millimeter activity for 13 days following its 
intial outburst. And whi le it achieved a si milar peak flux after 
distance correction (see Furuva et al. 2003, Fig. 2), the outburst 
activity does not appear to be periodic (like in the case of V773 
Tau A) since only a single radio flare was detected in the course 
of a 7-month monitoring campaign (lFellietal.1119931 As for 
ORBS, the radio flare shares a similar rise time a nd remains ac- 
tive s everal hours later when the track ended (see Forbrich et al. 
2008, Fig. 1). Unfortunately, little more is known about the 
deeply embedded YSO or its flare activity. Further comparison 
between these scenarios and DQ Tau is difficult, given that we 
are unable to constrain the emission polarization or its extent. 

Instead, from the millimeter perspective, the activity occurs 
within the window when overlapping magnetospheres from the 
binary components are expected to interact and produce flares 
like that observed on 19 April. And, in the case of V773 Tau 
A, similar star-star magnetospheric interactions were responsi- 
ble for a nearly identical 3 mm flare profile. However, pulsed 
accretion flows in the DQ Tau system also show a periodicity 
coinciding with periastron. And while magnetospheric interac- 
tion is likely to explain the millimeter flare, the magnetospheres 
do not possess the necessary energy requirements to power the 
optical brightenings with the frequency an d duration observed at 
optical wavelengths (M athieu et al.ll 1 997b . 

DQ Tau also possesses an additional component, absent in 
the V773 Tau A system, which must be considered: a circumbi- 
nary disk. Increased accretion at periastron inherently requires 
inclusion of a star-disk interaction. When stars accrete, material 
is thought to be swept from the (circumbinar y) disk, along th e 
magnetic field lines, and onto the protostars dShu et al.l [1994). 
These field lines can be stretched or compressed during accre- 
tion events, changes that can also produce flares. Furthermore, 
in accretion theory, the balance between gas pressure and mag- 
netic pressure often delineates the inner disk boundary. But 
when a stable accretion disk is disrupted (due to the dynami- 
cal motions of a binary) and not available to compress the star- 
disk field lines, radio flaring may occur more easil y, and over 
larger distances, as the field lines roam more freely (Bas ri et all 
1997). Increased outflow signatures during optical brightenings 
are consistent with increased field alignment, if magnetic fields 
drive the outflows, and if pulsed accret ion can organiz e the lines 
temporarily to power those outflows dBasri et al.ll 1997b . 

Therefore the success with each orbit of a dynamically- 
driven pulsed accretion mechanism could inhibit the incidence 
of millimeter activity (through compression of the field lines), 
resulting in an anti-correlation between optical and millimeter 
flaring. Alternatively, if the re-connection events provide a nec- 
essary pathway for accretion processes, then the haphazard na- 
ture of magneti c fields easily could explain the variability in the 
optical flaring (Ma thieu et al.|[T997l) and the occurrence of mil- 
limeter flares and optical brightenings is correlated. Finally, the 
optical brightenings and the millimeter flares may be indepen- 
dent, only sharing the common periodicity of the binary orbit. 

Simultaneous monitoring of this system at both optical and 
millimeter wavelengths is important to investigate whether mag- 
netospheric and dynamical interactions between the stars around 
periastron are strongly related. Ultimately, core fragmentation 
during star formation may favor multiple-star systems like DQ 
Tau. Binaries comprise 65% or more of the stellar popula- 
tion i n the middle of the main-sequence dDuquennov & Mavorl 
1 199 lb. And because plan ets have been detected around binaries 
(Egg enberger et alll2004l) . there is an established need to study 



these systems and the impact of magnetospheric interactions and 
variable accretion processes on disk evolution. 

6. Conclusions 

The PMS spectroscopic binary, DQ Tau, has a highly eccen- 
tric orbit that should allow the stellar magnetospheres of the 
two components to interact at each periastron passage. We 
serendipitously observed evidence for this activity in a flare 
at 3 mm, the first detected around a CTTS. The maximum in- 
tensity recorded (355 mJy beam -1 ) was 27 times the quiescent 
value (13.3 mJy beam -1 ). The observed peak occurred roughly 
7.6 (±18.2) hours before periastron at $ = 0.98+0.05, lasted 
several hours, and coincided with the average orbital phase 
(O = 0.9) of the well-documented, periodic optical brightenings. 
Whether a simultaneous optical flare occurred is unknown. 

The system's optical brightenings are consistent with the the- 
ory for pulsed accretion from the circumbinary disk onto the 
protostellar surfaces, regulated by the binary orbit, and perhaps 
assisted by magnetospheric interactions. The role of the magne- 
tospheres during accretion events can now be studied observa- 
tionally using methods established by Mas si et ail d2006l [2008) 
in the analysis of V773 Tau A, the first example of coronae in- 
teractions between protostars where the authors were able to re- 
solve the extent of the emitting region and the emission mech- 
anism using very long baseline interferometry and polarization 
measurements. A similar observing program aimed at DQ Tau 
offers a unique opportunity to study star-disk interactions and 
re-connection events directly, and to test if magnetospheric in- 
teractions are necessary to facilitate periodic accretion. 

Our findings also provide a small caution for millimeter flux 
points in SEDs that could contain unrecognized flare contribu- 
tions if measured only once, with the time consuming nature of 
observations at these wavelengths leading to a paucity of follow- 
up measurements. 
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